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Abstract 

Markovian evolving graphs are dynamic-graph models where the links among a fixed set of nodes 
change during time according to an arbitrary Markovian rule. They are extremely general and they 
can well describe important dynamic-network scenarios. 

We study the speed of information spreading in the stationary phase by analyzing the completion 
£\j | time of the flooding mechanism. We prove a general theorem that establishes an upper bound on 

flooding time in any stationary Markovian evolving graph in terms of its node-expansion properties. 

We apply our theorem in two natural and relevant cases of such dynamic graphs. Geometric 
Markovian evolving graphs where the Markovian behaviour is yielded by n mobile radio stations, 
with fixed transmission radius, that perform independent random walks over a square region of the 
(f) ' plane. Edge- Markovian evolving graphs where the probability of existence of any edge at time t 

depends on the existence (or not) of the same edge at time t — 1. 

In both cases, the obtained upper bounds hold with high probability and they are nearly tight. In 
fact, they turn out to be tight for a large range of the values of the input parameters. As for geometric 
Markovian evolving graphs, our result represents the first analytical upper bound for flooding time 
on a class of concrete mobile networks. 
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1 Introduction 



Markovian evolving graphs and Flooding. Graphs that evolve over time are currently a very 
hot topic in computer science. They arise from several areas such as mobile networks, networks of users 
exchanging e-mail or instant messages, citation networks and hyperlinks networks, peer-to-peer networks, 
social networks (who-trust-whom, who-talks-to-whom, etc.), and many other more [21 [TCI HU1 |21?1 |2~T1 [501 

ED. 

Markovian evolving graphs are a natural and very general class of models for evolving graphs intro- 
duced in [2]. In these models, the set of nodes is fixed and the edge set at time t stochastically depends 
on the edge set at time t — 1: so, we have an infinite sequence of graphs that is a Markov chain. It is 
important to observe that, on one hand, these models make the underlying mechanism of how the graph 
evolves explicit; on the other hand, they are very general since, by a suitable choice of the matrix transi- 
tion probability yielding the graph Markovian process, it is possible to model several important network 



scenarios such as faulty-networks and geometric-mobile networks (such scenarios will be described later). 

In 2] , the hitting time and cover time of random walks in some specific cases of Markovian evolving 
graphs have been analytically studied. We instead investigate the speed of information spreading on 
general Markovian evolving graphs. Reaching all nodes from a given source/initiator node is typically 
required to disseminate or retrieve information: this task is performed via suitable protocols that aim 
to achieve low delay and message overhead. However, when the network topology is highly dynamic 
and unknown, (e.g. unstructured peer-to-peer networks, faulty /mobile networks, etc), it is very hard to 
design efficient protocols and, as a result, the flooding mechanism is often adopted [5J [TCJ H7J [55] . In the 
flooding mechanism, any informed node (i.e. any node that has the source message) always sends the 
source message to all its neighbors. So, the source is informed since the beginning and, clearly, any other 
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node gets informed at time step t iff any of its neighbors (w.r.t. the edge set at time t) is informed at 
time step t — 1. 

The completion time of the flooding mechanism (shortly flooding time) is the first time step in which all 
nodes of the network are informed. 

It is important to observe that flooding time of a dynamic network may largely differ from its diameter: 
for instance, it is easy to construct an n-node mobile network over a finite square that has, at every time, 
diameter D = 3 while its flooding time is 0(n). In general, any diameter bound for a given dynamic 
network implies nothing about its flooding time but the fact that the latter is finite. Flooding time in 
fact represents the "natural" lower bound for broadcast protocols in dynamic networks. For this reason, 
flooding is often used in order to evaluate the relative efficiency of alternative protocols, especially in 
networks with unknown dynamic topology [5J 1161 129] . 

Our results. We study flooding time in stationary Markovian evolving graphs, i.e., when the initial 
graph is random with a stationary distribution of the underlying Markov chain pQ . In network mobility 
simulation, this corresponds to the important concept of perfect simulation (see |25l 17]). 

We prove an upper bound on flooding time in any stationary Markovian evolving graph. This upper 
bound is expressed in terms of the parameterized node-expansion properties satisfied by the stationary 
graphs. As far as we know, this is the first analytical result on the speed of information spreading in so 
general dynamic models. 

We then show the tightness (so the "goodness") of this bound in two relevant concrete network scenarios: 
geometric Markovian evolving graphs (in short, geometric- MEG) and edge- Markovian evolving graphs (in 
short, edge- MEG). 

Geometric Markovian evolving graphs. We consider a model of evolving graphs that is based on node 
mobility. It is the discrete version of the well-known random-walk model |19l l6l fT3l 122) . In this model, 
denoted here as geometric- MEG ', nodes (i.e. radio stations) move over a region of the plane (typically a 
square region) and each node performs, independently from the others, a sort of Brownian motion. At 
any time there is an edge (i.e. a bidirectional connection link) between two nodes if they are at distance 
at most R (typically R represents the transmission range). We make time discrete and consider a square 
grid of arbitrary resolution as a node support-space (see Section [3] for details) . This model can also be 
viewed as the walkers model [13] on the square grid. 

It is important to observe that geometric-MEG yield stochastic dependency among the dynamic edges, 
i.e., the probability of an edge depends on the existence of other edges. 

The impact of mobility in information spreading has been the subject of several papers over the last years. 
However, only few analytical results are currently available. In |18) . some bounds on the network capacity 
(i.e. the number of received packets) has been analyzed on a mobility model that is not explicitly defined. 
In |24j , the authors analyze the broadcast time over a restricted mobility model. In this restricted model, 
at every time step, the position of each node is selected independently at random inside a disk that is 
fixed at the starting time. Observe that, in this restriction, there is no stochastic dependence between two 
consecutive node positions: the model is significantly far from the random walk model. Then, the same 
work provides some experimental results for the random walk model. Finally, in [22], the speed of data 
communication between two nodes is studied over a class of Random- Direction models yielding uniform 
stationary node distributions (including the random walk model with reflection). They provide an upper 
bound on this speed that can be interpreted as a lower bound on flooding (routing) time when the mobile 
network is very sparse and disconnected (so, differently from our result, they consider geometric-MEG 
under the connectivity threshold). Their adopted technique based on Laplacian transform of independent 
journeys strongly departs from ours and it cannot be extended to provide any upper bound on flooding 
time. Further related analytical results that have been obtained after the conference version of our work 
are discussed in Section [3] 

We first prove that stationary geometric MEG, yielding connected graphs, satisfy certain parameterized 
node-expansion properties. We then apply our general result and achieve an upper bound on flooding 
time. The obtained bound is shown to be tight whenever flooding time is f2(loglogn). Informally 
speaking, this happens whenever (i) the transmission radius is not "almost" equal to the diameter of 
the square region and (ii) the maximal node speed is less than the message-transmission speed. Both 
assumptions are satisfied by most of real mobile networks. In general, our upper bound is thus at most 
an O (log log n) additive factor larger than the optimum. 

Further mobility models. The node-expansion properties of geometric-MEG are mainly due to the fact 
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that the stationary distribution of node positions is almost uniform. In this paper, we provide formal 
results and proofs only for flooding in geometric-MEG. However, our expansion technique can be applied 
to any mobility model yielding a uniform or almost uniform stationary distribution of node positions. 
Several variants of the random waypoint model, one of the most commonly used mobility models (231 
l6l 125). enjoy this uniformity property. Among the others, we mention the random-direction model with 
reflection (also called the billiard model) [31 [251 [25], the random waypoint on a torus [T51 |2T)1 [231 |2"B] 
and the random waypoint on a sphere [25] . Furthermore, the uniformity property is also satisfied by the 
walkers model on a toroidal grid [14] . 

To the best of our knowledge, our results are the first analytical bounds on flooding time for natural 
and concrete models of mobile networks. 

We finally remark that our flooding analysis does not take care of the interference problem in message 
transmissions: this is typically managed at the MAC layer of a wireless network architecture [3J [TU]. The 
impact of message interferences in geometric-MEG is a further interesting issue which is out of the scope 
of our work focussing instead on dynamic-topological properties of MEG. 

Edge-Markovian evolving graphs. In several network scenarios, there is a strong dependence between the 
existence (or the absence) of a link between two nodes at a given time step and the existence (or the 
absence) of the same link at the previous time step. Important examples of this behavior arise in faulty 
communication networks, peer-to-peer network^], and social networks. 

We thus consider edge-MEG, special Markovian evolving graphs, recently studied in [51 0], which 
are a time-discrete version of the reciprocity graph model introduced in the context of evolving social 
networks |32j . At every time step, every edge changes its state (existing or not) according to a two-state 
Markovian process with probabilities p(n) and q(n) where n is the number of nodes. If an edge exists 
at time t then at time t + 1 it dies with probability q(n) (i.e. death-rate). If instead the edge does not 
exist at time t, then it will come into existence at time t + 1 with probability p(n) (i.e. birth-rate). For 
brevity's sake, functions p(n) and q(n) will be simply denoted as p and q, respectively. Observe that 
setting q = 1 — p yields (time-independent) dynamic random graphs studied in |10j to model dynamic 
radio networks and in [5] to model epidemic biological processes; here links, at every time, are chosen 
independently at random. So, edge-MEG are (in turn) a wider and more realistic class of dynamic 
random graphs. Observe that when < p, q < 1, the stationary distribution is unique. 

Similarly to the case of geometric-MEG, we first prove that stationary edge-MEG, yielding connected 
graphs, satisfy certain parameterized node-expansion properties. Thanks to these properties, we can 
apply our general result and achieve an upper bound on flooding time. The obtained bound is shown to 
be tight whenever flooding time is fi(loglogn): this includes, for instance, the relevant case where the 
expected node-degree is O (poly log n). In general, our upper bound for edge-MEG is thus at most an 
O(loglogrt) additive factor larger than the optimum. 

In [5J, the maximal flooding time has been studied in edge-MEG with respect to any initial prob- 
ability distribution. In that paper, in fact, almost tight bounds for the worst-case flooding time have 
been derived. However, those results do not say whether flooding can be (significantly) faster in sta- 
tionary edge-MEG. Interestingly enough, our stationary bound implies that, whenever the birth-rate p 
is 0(l/n 1+c ) and the death-rate q is 0(np/ log n), there is an exponential gap between the stationary 
case and the worst-case. An exponential gap also holds whenever p — 0(log n/n) and q = 0(py/n) (for 
instance, set q — polylog n/n). 

Organization of the paper. In Section [21 we prove our upper bound for flooding time in general 
Markovian evolving graphs. The results for geometric-MEG and edge-MEG are described in Sections [3] 
and 21 respectively. Finally, further related analytical results (obtained after the conference version of 
our work) and some open questions are discussed in Section [S] 

1 Notice that, in some of these settings, there is an underlying physical network that supports the abstraction of point- 
to-point communication. 

2 Hence, any inequality p ^ (^)b(n) means that p(n) is eventually not larger (not smaller) than b(n). The same holds 
for q = q(n). 
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2 Markovian evolving graphs: the general theorem 

Through this paper, the set [n] = {1, . . . , n} will represent the set of n nodes. Let G = ([n], E) be a 
graph and I C [n] be a subset of nodes. We denote by N(I) the out-neighborhood of /, i.e. 



Given a source node s € [n], the flooding process can be represented by the sequence {I t C [n] : t e N} 
where 7 t is the subset of informed nodes defined recursively as follows 



Notice that the subset N(I t ) refers to the graph at time step t. Let T(s) be the first time step such that 
all nodes are informed. The flooding time is the maximum T(s) over all possible choices of source s. 

Definition 2.1 (Markovian evolving graph) Let G be a family of graphs with the same node set [n]. 
A Markovian evolving graph M. = {G t : t e N} is a Markov chain with state space G. 
A stationary Markovian evolving graph is a Markovian evolving graph M = {G t ■ t € N} such that Gq 
is random with a stationary distribution of M. 

The following definition concerns a sort of parameterized node-expansion. This is a key-ingredient, in 
our analysis of flooding in Markovian evolving graphs, to cope with the difficulties due to the stochastic 
dependence. 

Definition 2.2 (Expander) A graph G — ([n],E) is a (h, /c)-expandcr if, for every set of nodes I C [n] 
with \I\ ^ h, it holds that \N(I)\ > k\I\. 

The above definition naturally extends to random variables and their probability distributions. 

Definition 2.3 (Expander II) Let X be a random variable with values in a family of graphs with the 
same node set [n]. Then X is a (h, /c)-expander with probability p if 



In this case, we also say that the probability distribution of X yields an (h, k)- expander with probability 
P- 

We are now able to provide our main result for general stationary Markovian evolving graphs. We first 
show a lemma that connects the parameterized expansion of a (deterministic) evolving graph with its 
flooding time, then we use it to prove our theorem on the flooding time of stationary Markovian evolving 
graphs. In what follows, all logarithms are in base e. 

Lemma 2.4 (Flooding and Expansion: Deterministic case) Let Q = {G t ■ t e N} be an evolving 
graph (i.e. a sequence of graphs with the same node set [n]). Suppose an increasing sequence 1 = ho ^ 
hi < ■ ■ ■ < h s = n/2 and a non-increasing sequence ki ^ • • • ^ k s of positive real numbers exist such 
that, for every t £ N, graph Gt is a (hi, ki)-expander for every i = 1, . . . , s. Then the flooding time of Q 



Proof. Let m t be the number of informed nodes at time step t, at the beginning mo = 1. For i = 1, . . . , s 
let Tj be the first time step such that the number of informed nodes is larger than hi, 



If t is a time step such that < m t < hi for some i — l,...,s, then, since the graph G t is a 

(hi, fci)-expandcr, it holds that the number of informed nodes at the next time step is 



N(I) = {v e [n] \ I : {u, v} e E, for some u e 1} 



J t+ i = I t UN (I t ) 



P (X is a (h, k)-expander) > p 




Ti = min{t € N : m t > hi} 



m t+ i > (1 + ki)m t 
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Hence, after t' of such time steps it holds that 

m t+t < >{l + hf'mt > (1 + hfhi-i 
So the number of time steps between T,_i and Tj is at most 

log(/lj//l;_i) 



7i. — T,;_i 



log(l + kt) 



If W^T 1 = 0(1) then T t = O ( W&Y 1 ) ■ " '"f^ = °« thcn (1 + ^-i » 



iog(i+fei) — "v-v uncii -t, — ^ ^ iog(i+fei) /■ iog(i+fcj) 

that is in just one time step the number of informed nodes jumps from m t < /i, to m t+ i much larger 
than hi. Let j be the index such that hj < (1 + ^ /ij+i (if none of such index exists it means 

m t +i > n/2, so T s — Tj_i = 1), then it holds that 

iog(fej + i//ti_i) _ iog(/i£//i£_i) ^ iog(^/^_i) 

^ Iog(l + fci) log(l + fti) ^ log(l + fo) 

In the last inequality we used that ki ^ fcj for £ ^ i. Hence we can bound 



Tj — Tj_i < 



log(/lj + l/fej_l) 

log(l + fci) 



° E 



log(h e /he-i) 
log(l + fc*) 



Finally, by summing up the contributions of all the considered time intervals we have 

T„ = 



\ - bg(fti/fei-l) \ 

2, io g (i + fc 4 ) J 

Once we have at least n/2 informed nodes, then a symmetric argument holds. Indeed, consider the 
number fh t of non-informed nodes at time step t. Observe that the neighbors, in graph G t , of such 
nodes were not informed at the previous time step t — 1 (otherwise at time step t they would have 
informed their neighbors in graph Gt). Let i G {0, 1, . . . , s} be the index such that hi^\ < m t ^ hi, since 
G t is a (hi, h) -expander the number of such neighbors is at least kifh t , so the number of non-informed 
nodes at the previous time step t — 1 were at least 

fht-i > (1 + ki)m t 

In other words, the number of non-informed nodes follows the same growth rate, when there are at 
least n/2 informed nodes and we look at the system going backward in time, of the number of informed 
nodes, when there are less than n/2 informed nodes and the time moves forward. Hence, to go from n/2 

informed nodes to n informed nodes it takes further O log ( /l ' / ^'~ 1 - ) ) time steps. □ 



W=l log(l+fci) 

As usual, we say that an event (E)(n) holds with high probability (for short w.h.p.) if P (() (E)(n)) ^ 
1 - 1/n. 

Theorem 2.5 Let M = {G t : t e N} be a stationary Markovian evolving graph. Assume an increasing 
sequence 1 = ho ^ hi < ■ ■ ■ < h s = n/2 and a non-increasing sequence ki ^ • • • ^ k s (for any s < n/2) 
of positive real numbers exist such that, with probability 1 — l/n 2 ; for every i = l,...,s the stationary 
distribution of M. yields an (hi ,ki)- expander. Then the flooding time of M is w.h.p. 

\fy Iog(l + fci) ) 

Proof. For t = 0, 1, . . . define the event 

£t = { Gt is a (hi, /^-expander for every i = 1, . . . , s } 
By stationarity hypothesis we have that P (£<) ^ 1 — 1/n 2 for every t. Now consider the event 
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and observe that from Lemma |2~H it follows that HtLo ^* — hence T C |J"=o ^* an d we have that 

P(J) ^P ' 1 m " ^ ' 




□ 



An easy consequence of Theorem 12.51 is the following 

Corollary 2.6 Let M. = {Gt ■ t € N} be a stationary Markovian evolving graph. Assume a non- 
increasing sequence k\ ^ • • • ^ fc„/ 2 of positive real numbers exists such that, with probability 1 — , for 
every i = 1, . . ., n/2 the stationary distribution of M yields an (i, ki)-expander. Then the flooding time 
of M. is w.h.p. 

(n/2 

T.1 



—J i log(l + h) 

3 Geometric Markovian evolving graphs 

We introduce a model of dynamic graphs that is a discrete version of the random walk mobility model for 
radio networks [6]. In the latter model, nodes (i.e. radio stations) move on a bounded region of the plane 
(typically a square region) and each node performs, independently from the others, a sort of Brownian 
motion. At any time there is an edge (i.e. a bidirectional connection link) between two nodes if they 
are at distance at most R (typically R represents the transmission range). In our model we discretize 
time and space. We choose to keep constant the density (i.e. the ratio between the number of nodes 
and the area) as the number n of nodes grows. The node region is a square of side y/n and the density 
equals to 1 . This choice is only for the sake of simplicity and all the results can be scaled to any density 
S(n) (see Observation 13.31) . The nodes can assume positions whose coordinates are integer multiple of 
a sufficiently small resolution coefficient e > 0; in the sequel, we always assume that e ^ 1 and e < R. 
Formally, nodes move on the following set of points 

L n>e = {(ie,je) | i,j G N A i,j ^ ^-} 

e 

At any time step, a node can move to one of the positions of L n>e within distance r from the previous 
position. The positive real number r is a fixed parameter that we call move radius. It can be interpreted 
as the maximum velocity of a node[f|. Formally, we introduce the move graph M„ jr c = {L n ^, E n ^ rje ), 
where 

E n ,r,e = {(x,y) | x,y e L„, e d(x,y) < r} 

and d(-, ■) is the Euclidean distance. A node in position x, in one time step, can move in any position 
in F(x), where T(x) = {y | (x, y) £ £Vi.r, e }. The nodes are identified by the first n positive integers [n]. 
The time-evolution of the movement of a single node i is represented by a Markov chain {Pi t t ; t € N} 
where P^t are random variables whose state-space is L nj£ and 

In other words, P^ t is the position of node i at time t. Thus, the time-evolution of the movements 
of all the nodes is represented by a Markov chain V(n,r,e) = {Pt : i £ N} whose state-space is 
L n ^ t x x • • • x L„ j( (n times) and 

Pt — (Pl,t,P2,t, ■ ■ ■ , Pn,t) 

Let us fix a transmission radius R > 0. A geometric-MEG is a sequence of random variables Q (n, r, R, e) = 
{G t : tefi} such that G t = {[n],E t ) with 

E t = {{i,j)\d(P itU P jtt )^R} 

From a formal point of view, geometric-MEGs are not Markovian evolving graphs according to Defini- 
tion 13.11 In order to include them, we need a slight generalization of the definition 

3 Indeed, a node can run through a distance of at most r in a unit of time. 
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Definition 3.1 (Markovian Evolving Graph II) Let G be a family of graphs with the same node set 
[to]. A Markovian evolving graph Q = {Gt : t £ N} is a sequence of random variables with state space G 
and such that there exist both a Markov chain X — {X t : t £ N} and a function f so that Gt = f{X t ). 
A stationary Markovian evolving graph is a Markovian evolving graph Q = {Gt '■ t £ N} such that Gq is 
random with a stationary distribution of X translated by f. 

Theorem 12.51 easily extends to the above generalized definition by straightforward arguments. 

As for the stationary case, standard results of Markov chain theory (see pQ) easily imply that the (unique) 

stationary distribution ^ of Markov chain {Pi,t t £ N} is 



Ti; (x) 



£ ye w l r (y)l 



Notice that 7Tj is almost uniform since, for any two positions x and y, the values 7Tj(x) and 71; (y) can differ 
by at most a constant factor. Moreover, the stationary distribution of V(n, r,e) is the product of the 
independent distributions 7T; for all i £ [to]. We say that a geometric- MEG G{n,r,R,e) = {Gt : t £ N} 
is a stationary geometric-MEG if the underlying Pq is random with the stationary distribution of the 
Markov chain V{n,r,e) = {P t : t £ N}. Notice that if Q(n, r, R, e) = {G t : t £ N} is a stationary 
geometric-MEG then all random variables Gt are random with the same probability distribution that 
we call stationary distribution of Q(n, r, R, e). 

Stationary geometric-MEG enjoy of the following expansion properties. 

Theorem 3.2 Ifr^O and cyTogro ^ R ^ \fn for a sufficiently large constant c, then constants a, fj > 
exist such that, with probability 1 — -4-, the stationary distribution of Q(n,r, R,e) yields: 

• A (h, a ht)- expander for 1 ^ h ^ aR 2 ; 

• A (h, (3^=) -expander for aR 2 < h ^ to/2. 

Observation 3.3 For general node density 5(n), Theorem ] 3. OA holds under the scaled assumption R 
cy/logn/8(n). 

Proof. Let m = \y/En/K\ . Consider the partition of the square y/n x y/n into m x m congruent sub- 
squares, called cells. Every cell can be identified by the pair of indices (i,j), for 1 ^ i,j ^ m, such that 
i is the index of row and j is the index of column of the cell. Let Cij be the subset of the points of L n fL 
that fall into the cell Notice that the side length I of a cell satisfies 

R/(VE+l) R/Vb 

Thus, any point of a cell is at distance less than R from any point of a side-by-side adjacent cell. 
Through the following, we assume that the positions of the nodes are random with the stationary 
distribution of the Markov chain V(n, r, e). Moreover, we say that a node belongs to a cell whenever its 
position belongs to the cell. Let N^j be the random variable counting the number of nodes in cell ctj. 
Now, we prove a simple but crucial claim. 



Claim 1 If e 1 and R cy/Togn for a sufficiently large constant c, then a constant X 1 
exists such that, with probability 1 — it holds that, for every 1 ^ i, j ^ m, 

^- < Nij < XR 2 

Proof. Firstly, consider a fixed cell (i,j). For every u £ [to], let X u be the {0,1} random 
variable that is 1 iff node u is in the cell Cjj. Clearly, these are independent random variables 
and it holds that N^j — YlueM ^ s ^ or probability distribution of X u , we have that 

P(X» = 1) = 
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Since 

n (x) = irWl 
4 j ~ E, eL „. £ |r(y)l 

it is easy to see that, if e is sufficiently small (say e ^ 1) then there is a constant 7^1 such 
that, for every x € I„ je , it holds that 

1 7 

^ 7T u (x) ^ 



l\L n ,e\ \Ln,e\ 

This implies that 

-M- <p(x„ = i) < iM 

By taking into account the side length of the cells, it is easy to verify that 

lOn \L n ,e\ 5n 

It follows that 

< P(X U = 1) < ^ and ^[N i3 ] < ^ 
107n v ' 5n IO7 1 Jj 5 



Now, in virtue of the ChernofFs bound [37], if i? ^ c-^/log n, for a sufficiently large constant 
c, then a constant A ^ 1 exists such that 

A 

with probability at least 1 — ^. Since the number of cells is less than n, a simple application 
of the union bound proves the thesis of the claim. □ 



Let B be the event that occurs when, for every 1 ^ i, j ^ m, 

sC N U1 \R 2 



If 
A 



where A is the constant of Claim [T] We now prove event B implies the expansion properties stated in 
the thesis of the theorem. 

Claim 2 // event B holds then the graph induced by R and by the positions of the nodes is a 
(h, a^-)-expander for 1 ^ h ^ aR 2 , where a — 1/(2A). 

Proof. Let I C [n] be such that |/| ^ aR 2 . Consider a node u in I and let Cij be the cell that 

contains u. Since B holds, Ni.j ^ ^j-. All the nodes in Cij are adjacent to u. Thus, there are 
at least Nij — \I\ nodes that are adjacent to u and that are not in 7. It follows that 

\N(I)\ > Ntj - \I\ > ^ -aR 2 > ^ = aR 2 

In other terms, \N{I)\ ^ af^\I\. □ 

Claim 3 // event B holds then the graph induced by R ^ y/n and by the positions of the 
nodes is a (h, (3-^) -expander for aR 2 ^ h ^ n/2, where (3 = 

Proof. Let I C [n] be any subset of nodes with |/| ^ n/2. We say that a cell is black if it 
contains at least one node in I. Let B be the random variable that is the set of black cells. 
Let J be the random variable defined as follows 

J = {u <G [n] j u $l I A 3c € B : node u belongs to c} 
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Now, two cases are possible: either |J| f3Ry/\I\ or not. Firstly, suppose that | J| ^ f3R\f\I\. 
Since every node in J is in a black cell, it holds that J C N(I), and thus 

iiv(/)i > iji > 

In other terms, iV(7) ^ /3— ^=|/| and the expansion property is proved. 



Consider now the case \J\ < (3Ry/\I\. We say that a row (column) of cells is black if all the 
cells of the row (column) are black. Similarly, we say that a row (column) is white if all the 
cells of the row (column) are white. A row (column) that is neither black nor white is said to 
be gray. Notice that a gray row (column) contains at least two adjacent cells such that one 
is non-black and the other is black. Let B r and B c be, respectively, the number of black rows 
and the number of black columns. Three cases may arise. 

[B r ^ 1]: Observe that in this case all the columns are either black or gray. Let Y be the 
number of gray columns. It holds that 

\B\ 

Y^m-B c ^m- 1 —^- 
m 

Since event B holds, the number of nodes in non-black cells is bounded by XR 2 (m 2 — \B\) 
and thus 

XR 2 (m 2 - \B\) > n-\I\-\J\ > n- \I\ - /3R\/\T\ 



It follows that 



LB I sS m 2 - 



2 n - \I\ - PRy/\T\ 



XR 2 

By combining this bound with the previous bound on Y we obtain 



Y > n-\I\-/3R^/\I\ > n-\I\-mV\I\ 



XR 2 m X2V$n~R 
where the last inequality follows from m — \\/5n/ R~\ and R ^ y/n. 

Observe that every gray column contains at least one non-black cell that is adjacent to a 
black cell. So, all the nodes belonging to those non-black cells are included in N(T). Since 
event B holds, it follows that 

\N(I)\ ^Y^-^R 

Now, recalling that /? = g^-, |/| < n/2, and R ^ y/n, it is easy to verify that 




X 2 2V5n 

It follows that \N(I)\ > j3R\f\T\ and the expansion property holds. 
[B c ^ 1 (and B r = 0)]: This case is symmetric to the previous one. 

[B r = and B c = 0]: In this case, all the rows and columns are either gray or white. Let Y r 
and Y c be the number of gray rows and the number of gray columns, respectively. Since there 
are neither black rows nor black columns, it must be the case that every black cell belongs 
to both a gray row and a gray column. As a consequence it holds that Y r ■ Y c ^ \B\. 
Without loss of generality, assume that Y r ^ Y c . Then Y 2 ^ \B\ and thus Y r ^ Since 
event B holds and every gray row contains a non-black cell adjacent to a black one, it holds 
that 

\N(I)\ > Y r ^ > ,/\B\^ 
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By using again the fact (implied by event B) that every cell contains at most \R 2 nodes, we 



have that \B\XR 2 ^ |/| and thus y/\B\ > ^H. It follows that 



\N(I)\ > Ml > PRy/\T\ 

and the expansion property holds. □ 

Since, by Claim [TJ event B occurs with probability at least 1 i?, Claims [2] and [3] imply that also the 

expansion properties will hold with probability 1 ^. □ 

Thanks to the general bound given by Corollarv l2.61 the above expansion properties can be exploited in 
order to bound the flooding time in stationary geometric-MEG. 

Theorem 3.4 Let Q[n,r, R, e) be a stationary geometric-MEG. If r ^ and cy / Iog _ n ^ R ^ y^n for a 
sufficiently large constant c, then the flooding time of Q(n,r, R, e) is w.h.p. 

0^+loglogii 

Proof. From Theorem 13.21 the stationary geometric-MEG Q(n,r, R, e) enjoys, with probability 1 — 
the following expansion properties: 

• (h, a -^-) -expander for 1 ^ h ^ aR 2 

• /3 -^-expander for aR 2 ^ h ^ n/2. 

Thus, by applying Corollary |2.6l we obtain that flooding time is w.h.p. 

(aR 2 l n/2 ^ 

° [ £ fclog(l + a^) + hlo g (l + f3^) 

We now evaluate the above two sums separately. For the sake of convenience, set T = aR 2 . It holds that 

T T 

^Mog(l + f) " ^Mog(l + ^)(T + M 
This holds since (y+h) ^1/2 for h ^T. Moreover, 

T T 

<A 1 T T <A 1 T 

^ ih \og(l + I) (T + h) ~ (T + 1) log(T + 1) + ^ + |) (T + h) 

,T y 

< 1 + / — ; ; Tfr-dx 

h x{T + x)\og{l + l) 

= 1 + [- loglog(l + = loglog(T) + c 
where c is a constant. Therefore we have shown that 

aR 



E ,, 77~T flST - ©(loglogi?) 



Now consider the second sum. By using the inequality log(l + x) ^ we have that 
^ 1 ^ y/fr + lg < 1 + & ^ J_ 
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where the last inequality comes from inequality 



for h aR 2 . Moreover, it holds that 

1/2 



y -7= < / -7= < 



By combining the above inequalities we obtain 



n/2 1 + J- 

V - < 2 ^ 



that is, 

n/2 



□ 

We remark that the proof of the expansion properties of Theorem 13.21 only relies on the fact that the 
stationary distribution of node positions is almost uniform. In fact we can get the same expansion 
properties for any mobility model yielding a stationary distribution of node position that is uniform or 
almost uniform. As mentioned in the Introduction, several relevant mobility models enjoy this uniformity 
property. So, thanks to our Theorem l2.5l we can get an upper bound on flooding time similar to that of 
Theorem [331 

Next theorem shows a lower bound on flooding time in stationary geometric-MEG. 

Theorem 3.5 Let Q(n,r, R,e) be a stationary geometric-MEG. If r 0, then its flooding time is w.h.p. 



R + r 

Proof. Since the geometric-MEG is stationary, it is not hard to see that, w.h.p., at time there exist 
at least two nodes u and v that are at distance greater than y/n/2. Consider the flooding process with 
source node v. Let xo be the position of v at time 0. For any t, let dt be the minimum distance from 
Xo that node u has ever reached during the first t time steps. It is immediate to see that dt+i ^ d t — r. 
Since do ^ i/n/2, it holds that d t ^ y/n/2 — r ■ t. 

Let D t be the maximal distance from xo that any informed node has ever reached during the first t time 

steps. It is easy to see that D t +\ ^ D t + R + r. Since D = 0, it holds that D t < (R + r)t. 

Let t be the time step in which node u gets informed. It must be the case that D T ^ d T . It follows that 

{R + r)r D T > d T ^ y/n/2-r-T. 

It follows that r ^ ■y r n/(2(R + 2r)). Therefore, the flooding cannot be completed in less than (^^j 

time steps. □ 
By comparing Theorem 13.41 and Theorem 13.51 we obtain the following 



Corollary 3.6 LetQ(n, r, R, e) be a stationary geometric-MEG. Ifr = 0(R), and cyTog n ^ R ^ i og i " g n 
for a sufficiently large constant c, then the flooding time of Q(n,r, R, e) is w.h.p. 



e 



~R 



Under the very reasonable conditions of the above corollary, the general bound on flooding time in 
Markovian evolving graphs thus turns out to be asymptotically tight for stationary geometric-MEG. 
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4 Edge-Markovian evolving graphs 

We recall the model introduced in An edge- MEG M(n,p,q) = {G t : t £ N} is a Markov chain 

such that Gt — ([n],E t ) with 

where {Xt(e) : e S ('2')} are independent Markov chains with transition matrix 





( 


1 


M = 





1 -p p 




V 1 


q l-q 



Remind that p is the birth-rate and q is the death-rate and notice that an edge-MEG is a Markovian 
evolving graph according to Definition l2.il Observe that if < p, q, < 1 the Markov chains {X t (e) : t £ 
N} are irreducible and aperiodic; so there is a unique stationary distribution 



P + 1 P + 

Hence, the stationary distribution of M(n,p, q) is G„,p (i.e. Erdds-Renyi distribution in which each 
possible edge occurs independently with probability p) where here and in the sequel 

P 

p + q 

Stationary edge-MEG enjoy the following node-expansion properties. 

Theorem 4.1 Let Ai(n,p,q) be an edge-MEG such that p ^ c^f^ 1 for a sufficiently large constant c. 
Then, the stationary distribution of A4(n,p, q) yields, with probability at least 1— a (h, -expander 
for 1 ^ h \ and a [h, -expander for \ ^ h ^ ^. 

The proof of the theorem is a simple consequence of the expansion properties of the G„ iP model (see for 
instance [H]). We here provide a detailed proof for our specific setting. 

Lemma 4.2 Let p ^ c log " for a sufficiently large constant c. With probability 1 \ for G n ^p it holds 

that for any L C [n] with \L\ ^ ^, 

\I\np n 



\N(I)\ ^ min 



c c 



Proof. We first consider the case when \I\ ^ A and prove that, with probability at least 1 — 4^, it holds 



\N(L)\ > Then we consider the case 4 ^ |J| ^ || and prove that, with probability at least 1 — 4y 

it holds \N(L)\ 2. 
Let to = \I\ < A. F< 

X u = otherwise. Since P (X v — 1) ^ we have 



Let to = |J| ^ i. For any u E [n]\I consider the random variable X u so that X u = 1 if u € AT (7) and 



E[|iV(/)|]= ]T E[X tt ] = (n-m)m^> 

«e[n]\/ ^ P 



A 1 



From ChernofFs bound we get 
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Therefore 



P (31 C [„], 1 < |/| < l/p : \N(I)\ < < E P (l^WI < ^l 7 ^) 

/ C [n] 



1 < |/| < l/p 

LVpJ , x 



m— 1 x 7 

LVpJ 1 

E„ _£=A rrt 1 c-S c-12 

n n 4 < -n 4 < n 4 

m=l y 



And the thesis follows if we choose c ^ 20. 

Now consider the case where A ^ |/| = m ^ Notice that |iV(J)| ^ 7 if and only if there exists a set 
AC[n]\(JU N(I)) such that \A\ > n - m - 2. Hence 



P { 3IC [ n],\I\ = m:\N(I)\^-) = ^ m ^ ^ j(l-p) 
From the following inequalities 

• (W = < ("f m ) < (?) < M ? = e?1 ° S(eC) < e" 1 "^^ 1 ^) 

• (1 - p) m \ A \ ^ e _mp|j41 < e~ mj3 (™~^~ = ) = e - " 1 ™*^ - ^ 
by choosing c sufficiently large, we get 

P (3/ C [ n ], |/| = m : |JV(2)| < ^) ^ e^"^* - * -3 ^) sc: e ~* 

Hence 

Pf3JC[»],-<|J|<^:|iV(/)|<^) < £ p(3IC[n],|J|=m:|Ar(7)|<^ 

V P C/ m=Ll/pJ 

fn/21 

m=Li/pJ 

where the last inequality holds for sufficiently large n. □ 
The expansion properties of stationary edge-MEG, stated in Theorem l4.ll allow us to apply Corollarv l2.6l 
and thus get the following 

Theorem 4.3 Let A4(n,p,q) be a stationary edge-MEG such that p c ^f^ f° r a sufficiently large 
constant c. Then flooding time in M(n,p, q) is w.h.p. 

(log n 
— 2— + loglog(np) 
log(np) 

Proof. Thanks to Theorem 14. 1[ we can apply Corollary 12.61 with sequence 



22 for 1 «; i «; 

C 

^ for 



Thus we obtain that the order of flooding time is w.h.p. bounded by 

LVpJ -, r»/cl-i , n/2 



^ iloefl + ^) + ^ ilog(l + ^-) + E ilog(l + ^) 

i=l ^sV 1 ^ c ; i=|_l/pj + l 6V cjV i=[ n /cl sv 
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We now evaluate the above sums separately. For the first sum, by using Y^iLi \ ^ logm + 1, we have 

^»log(l + 2fi) log(l + f) \\o%k4)) 
For the second sum, by using log(l + x) logo;, we have 

r»/ci-i 1 r«/ci-i ^ „|- n / c ]_x 

^ ^ / rfa; 

i lnoYl -I- — 'I ^-^ 7'lncrli /,,,., -nine -2- 



,= Li/pJ+i ll0g(1 + ^ ) i=LVftl+i il0g H Av#J ^ 



rn/cl-l 

LVpJ 

For the third sum, we apply log(l + x) ^ + a;) for x < 1 and get 



r 7i I r™/ c l _1 

= -loglog— = C(loglog(np)) 

L cxUi/si 



n/2 n/2 , „ s «/2 n/2 

E E ^t«— = E ; + < E S + S 

= Tn/cl 6V Cl/ i=[n/cl Ci i=[n/c] V 7 i=r n /cl 



□ 



Next theorem gives a lower bound on flooding time in stationary edge- MEG. 



Theorem 4.4 Let A4(n,p,q) be a stationary edge-MEG such that p c^Tp 1 / or a sufficiently large 
constant c. Then the flooding time of M(n,p, q) is w.h.p. 



n 



logn 
log(np) 



Proof. Let A t be the random variable indicating the maximal node degree of Gt- Since the marginal dis- 
tribution of Gt is Erdos-Renyi G n ,p then, for a sufficiently large c (say c = 4), it holds that P (At > 2np) < 
1/n 2 . By applying the union bound the probability that a time step t < n exists such that A t > 2np is 
at most 1/n. Thus, the number of informed nodes at time step t < n is at most (2np) w.h.p., and this 
number is less than n/2 for t < . □ 

By comparing the upper bound of Theorem 14.31 and the lower bound of Theorem 14.41 we obtain the 
following 



Corollary 4.5 Let M{n,p,q) be a stationary edge-MEG such that c^^ 1 < p ^ "'° s ' og " , for a suffi- 
ciently large constant c. Then flooding time in A4(n,p,q) is w.h.p. 



e 



Jog_n_\ 
log(np) / 



5 Conclusions 

We showed that in geometric-MEG, under some conditions on the maximal node speed and transmission 
radius, node mobility has an almost negligible impact on flooding time: the latter turns out to be 
equivalent to the diameter of the static stationary graph. A similar fact holds for edge-MEG as well. 
After the conference version of this paper, an improved, dynamic version of our expansion technique has 
been derived in [TT] in order to obtain almost tight bounds for the flooding time of highly-sparse and 
disconnected geometric-MEG when the maximal node speed is larger than the transmission radius. So, 
in this case, mobility significantly speeds-up flooding time with respect to the static case. 

An important open issue is to provide analytical bounds for the flooding time of evolving graphs that 
are somewhat non homogeneous. Interesting examples arc the evolving graphs yielded by node performing 
random walks over highly-irregular support graphs and those yielded by nodes moving according the 
random-waypoint model over a non-convex, irregular region. 
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